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Abstract 

HFOLD (Higgs Full One Loop Decays) is a Fortran program package for 
calculating all MSSM Higgs two-body decay widths and the corresponding 
branching ratios at full one-loop level. The package is done in the SUSY 
Parameter Analysis convention and supports the SUSY Les Houches Accord 
input and output format. 
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A future high-energy e~^e~ hnear colhder will be the best environment for the pre- 
cise measurements of masses, cross sections, branching ratios, etc.. Experimental 
accuracies are expected at the per-cent down to the per-mille level. These must be 
matched from the theoretical side. Therefore higher order calculations are manda- 
tory. 

Solution method: 

This program package calculates all MSSM Higgs two-body decay widths and the 
corresponding branching ratios at full one-loop level. The renormalization is done 
in the DR scheme following the SUSY Parameter Analysis convention. The pro- 
gram supports the SUSY Les Houches Accord input and output format. 
Restrictions: 

Unusual features: 

Additional comments: 

Running time: 



1. Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is the most ex- 
tensively studied extension of the Standard Model (SM) of elementary par- 
ticles. Supersymmetry (SUSY) provides a solution to the so called hierarchy 
problem and furthermore, in the context of this work, it is a renormalizable 
theory. If the MSSM is realized in nature, supersymmetric particles will be 
produced at the LHC. However, even if SUSY is discovered, it will still be a 
long way to determine the parameters of the underlying model, which would 
shed light on the mechanism of SUSY breaking. A future high-energy e"'"e~ 
linear collider will be the best environment for the precise measurements of 
masses, cross sections, branching ratios, etc.. Experimental accuracies are 
expected at the per-cent down to the per-mille level [H, [i], l^l . These must be 
matched from the theoretical side. Therefore higher order calculations are 
mandatory. 

For the decays of the MSSM Higgs bosons, the one-loop corrections due 
to gluon and gluino exchange (SQCD) are known analytically, see e.g. 0, Is], 
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3, B , l ol • Full one-loop calculations were done e.g. in 0, M, ini, H Q E 



16l . llTl . Il8l . Il9| . For calculating the full (including electroweak corrections) 
one-loop decay widths automatic tools for generating all Feynman graphs, 
and subsequently the squared matrix elements, are strongly needed. 

There are a few program packages available for the automatic computa- 
tion of amplitudes at full one-loop level in the MSSM: FeynArts/FormCalc 



29|, SloopS [20|,|21| and GRACE/SUSY-loop j22|. SloopS and GRACE/SUSY- 
loop also perform renormalization at one-loop level. However, so far there 
is no publicly available code for the two-body Higgs decays at full one-loop 
level in the MSSM. Therefore, we have developed the Fortran code HFOLD 



23(1 . It follows the renormalization prescription of the SUSY Parameter Anal- 



ysis project (SPA) |25|] and supports the SUSY Les Houches Accord (SLHA) 



input and output format |2J]. The package HFOLD (Higgs Full One-Loop 
Decays) computes all two-body decay widths and the corresponding branch- 
ing ratios of the three neutral and charged Higgs bosons at full one-loop 
level. 

This paper is organized in the following way: First we shortly recapitulate 
the Higgs sector in the MSSM. Then we will discuss the renormalization used 
in the program. We will compare the total and partial decay widths of the 
Higgs bosons at the SPSla' point with existing programs. The last section 
will be the program manual. 



2. MSSM Higgs sector at tree-level 

2.1. Masses and mixing angles 

In the MSSM two chiral Higgs superfields with opposite hypercharge are 
necessary to keep the theory anomaly free. Two Higgs doublets are also 
necessary in order to give separately masses to down-type fermions and up- 
type quarks. 

The scalar components of the two complex isospin Higgs doublets 




represent eight scalar degrees of freedom (d.o.f.) and have hypercharges 
Y{Hi^2) = =F1- After spontaneous electroweak symmetry breaking, their neu- 
tral components receive vacuum expectation values, = vi and (if^) = 
V2- The absolute value v"^ = vf + f| can be determined from the measure- 
ments of e.g. and the SU(2) coupling g, but tan (5 = — remains a free 
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parameter. There remain five physical Higgs bosons, two neutral CP even 
ones, and and one neutral CP odd field and two charged Higgs 
bosons H^. The physical states and are mixtures described by the 
mixing angle a. The remaining three d.o.f. are 'eaten' by the longitudinal 
components of the now massive vector bosons Z° and W^. 

At tree-level only two free parameters describe the Higgs sector. In the 
MSSM usually the parameters m^o and tan/3 are chosen. The other three 
Higgs boson masses and the mixing angle a can be expressed at tree-level 
by rriz and mw, e.g. m'jf+ = m^r + m\o- Contrary to the SM, the Higgs 
self-interactions are completely fixed by EW parameters. At tree-level the 
mass of the lightest Higgs boson cannot be larger than mz- This value 
is already ruled out by LEP2. Fortunately, radiative corrections push the 
theoretical limit up to mhO<135 GeV with the leading contributions from 
top and stop loops proportional to m^/rriyy. 
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a) 

scalar 



c) 

scalar 



Figure 1: Four possibilities of two-body decays of a scalar particle 

2.2. Decay patterns and some properties 

As fermion number is conserved we only have four possibilities of Feynman 
graphs (at any loop level) for a two-body decay of a scalar: the decay into 
two scalars, into two fermions, into a scalar and a vector particle, and into 
two vector particles, see Fig. [H In the case of Higgs bosons the following 
decays are calculated: 
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77, gg, 'jZ^ (loop induced) 



(f) = h^, H^, and f = i'u e, /x, t, m, (i, c, s, t, 6, /' denotes the isospin 
partner to f, e.g. f = t, f = b, f and / denote the SUSY partners of 
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/ and /', x'^ and are the neutralinos and charginos, respectively. The 
Higgs bosons couple to fermions via their Yukawa couplings. Therefore, the 
branching ratio (BR) into top quark(s) is large, if the decay is kinematically 
allowed. The BRs oi hP ^ hh and to r+r~ are dominant, especially for 
large tan /3. The decays into the third generation sfermions may become 
dominant when they are kinematically possible. The decays into quarks and 
squarks can have large one-loop SQCD corrections. The decays into charginos 
and/or neutralinos can have significant one-loop contributions from the third 
generation (s)fermions depending on the gaugino/higgsino mixing. 

Decoupling limit: In case of m^o ^ mzo the masses of H^, A^, and 
become degenerate. 

This limit is already reached to a good approximation for m^o ~ 300 GeV. 
Furthermore, the {h^,H^) mixing angle can be expressed by a — /3 — 7r/2. 
Thus, the properties of the lightest Higgs boson are almost indistinguish- 
able from those of the SM Higgs boson. As a consequence, the couplings to 
the heavier Higgs bosons vanish at tree-level, e.g. the H^W~h^ coupling is 
oc cos(/3 — a) — 7- 0. 



3. Calculation at full one-loop level 

The definition of the MSSM parameters is not unique beyond the leading 
order and depends on the renormalization scheme. Therefore, a well-defined 
theoretical framework was proposed within the so-called SPA (SUSY Param- 



eter Analysis) project [25[. The "SPA convention" provides a clear base for 
calculating masses, mixing angles, decay widths and production processes. It 
also provides a clear definition of the fundamental parameters using the DR 
(dimensional reduction) renormalization scheme. These fundamental param- 
eters can be extracted from future collider data. The formulae for the wave 
function and mass counterterms (CTs) for sfermions, fermions and vector 
bosons in the on-shell scheme derived from their renormalization conditions 



can be found e.g. in [26|, l27|, |28 



The code of HFOLD is derived in the SPA convention in the general linear 
gauge for the IV^ and Z°-boson. All amplitudes are generated by using 
the tool FeynArts (FA) and the Fortran code is produced with the help of 
FormCalc (FC). For that purpose we imported all necessary formulae for the 
CTs into a FA model file. 
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The renormalized one-loop amplitude is the sum of the tree-level ampli- 
tude and the one-loop contributions, see Fig [21 

The tree- level couplings are given at the scale Q, implying that there are 
no coupling CTs. The DR scheme is defined by setting the UV divergence 
parameter A = 0. We however work with A 7^ and take for the coupling 
CTs only the parts oc A. In case the renormalized amplitudes are finite it is 
a proof for RGE invariance of the ordinary DR scheme. 




J^l _y\^Vertex J^CT _A^WFR 



Figure 2: One- loop renormalization procedure of a 1 to 2 process schematically 

The vertex corrections and all selfenergy contributions except the diago- 
nal wave function corrections can be directly calculated with FA/FC. 

Since there are many decay channels it was worthwhile to develop an 
automatic code generator at Mathematica level. First of all, it was necessary 
to work out all counterterms (in Mathematica form) for the whole MSSM. 
The idea is, not to have all MSSM couplings (which are more than 300 ones) 
at one-loop level hard coded in the MSSM model file of FA, but to calculate 
locally the amplitudes with the wave function and the coupling CTs (see 
Fig.E]). 

For each external particle we get a contribution to the wave function 
CTs amplitude by multiplying the bare fields with the corresponding wave 
function renormalization constants. The amplitude for the coupling CTs is 
obtained in the following way: First we calculate the tree-level amplitude, 
then we shift all tree- level couplings by their corresponding counterterms 6gi, 
Qi ^ Oi + Sgi and then take into account only terms linear in Sgi. 

The total two-body Higgs decay width can be written in one-loop approx- 
imation as 

Nc X kin X (lA^oP + 2 Re{MlMi)) , 
K,{mQ, mf, 777,2) 



r = 

kin = 
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with the totally symmetric Kdllen function k{x, y,z) = ^ {x ~ y — zY — Ayz 
and the color factor Nc = 3 for decays into quarks and squarks and Nc = 1 
for decays into other particles, respectively. 

M.I denotes the UV finite one-loop amphtude. The prefactor kin is a 
function of the on-shell masses of the incoming Higgs boson and outgoing 
particles only. Massless particles in loops can cause so-called infrared (IR) 
divergences in F. For this purpose, a regulator mass A for the photon and 
gluon is introduced. Adding then real photon or gluon radiation cancels these 
divergences. 



4. Input parameters 

HFOLD is designed to be applied to SUSY models like mSUGRA, GMSB 
or AMSB, where the low energy model parameters are given at some scale 
Q. The low energy spectrum is derived from a few parameters defined at 
a high scale using renormalization group equations. At the program start 
HFOLD reads the spectrum, where the Yukawa couplings, the gauge cou- 
plings gi,g2,g3, the soft breaking terms, the VEVs, m^o,tan/3, /i and the 
on-shell Higgs masses are taken as input parameters. The input parameters 
are understood as running parameters in the DR scheme at the scale Q. In 
loops we are free to use DR masses because the difference is of higher order 
in perturbation theory. Since our renomalization is done in the DR scheme 
the coupling counterterms contain only UV-divergent parts. Therefore we do 
not fix 6mw with Gpermi as input parameter. In the Higgs sector we use m^o 
and the running tan /3 as inputs. We can then simply derive the DR running 
Higgs mixing angle a at the scale Q. We do not take aeff as input parameter 
because we consider our calculation a self-consistent one-loop expansion. 



5. Resummation of tan/3 

The down-type fermions couple to the up-type Higgs doublet with radia- 
tive corrections by 

- y,H^ibb-y, A, cot PH^Jb. (1) 

The selfenergy A;, is proportional to tan /3 and can be enhanced for large val- 
ues of tan p. This term can be resummed (in the effective potential approach) 
by replacing the bottom Yukawa coupling jsoj with 
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The resummation can also be performed in the diagrammatic approach [31 . 
Different renormahzation schemes correspond to different choices of countert- 
erms. Therefore the analytic form of the tan /3 enhanced corrections depend 
on the chosen renormahzation scheme. In the on-shell scheme one takes the 
measured bottom mass as input parameter. The choice of 6mb fixes 6yb by 

nib r Smb 
yb = > (>yb = . (3) 

Vd Vd 

The quark mass counterterm 6mb is a source of tan /3-enhanced corrections. 
The selfenergy T,^^ {rnb) contains terms proportional to yb sin {3 and is there- 
fore tan /3 enhanced, 

E^^ = mbAb, (4) 
A, = Af + Af + Af . (5) 

In leading order this means : dnib = — = — mbe^tan /3. We write the bare 
Yukawa couplings as y^^^ = yb + ^Vt, where yb is the renormalized coupling 
and 6yb is the counterterm. The choice of 6mb fixes 6yb through 

Syb = = -ybebtan(3. (6) 

Vd 

The supersymmetric loop effects encoded in eb enter physical observables 
only through 6yb- Choosing e.g. a minimal subtraction scheme like the DR 
scheme for 6mb removes the tan /3-enhanced terms and there is nothing to 
resum anymore. Since we do not use the measured bottom mass as input, the 
resummation of tan /3 is absent in our approach. However, the resummation 
eq. ([2]) is implemented in the code and can be turned on. 

5.1. Gauge used 

The gauge fixing Lagrangian in the general linear gauge is given by 

with F+ = d^W^ + i^wmwG+, = ^^Z^' + izmzG'^.F^ = 7^^^, and 
F9 = -f^G^f". 
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The Higgs-ghost propagators are — C,v^v) and the vector-boson prop- 
agator reads 



Dv 2 -2 • 

— rriy 

The ^-dependent part is a product of two propagators leading to a (n+1)- 
point loop integral. Performing a decomposition into partial fractions, it can 
be split into a form with single propagators only, 

_ -^9'"' ^ ^ ( q'^q" q^q" \ 



q^ — my rriy \q^ — my q^ — ^rriyj 

We have implemented this second form into FA in order to check the gauge 
independence for W and Z. For the massless particles 7 and gluon we get 
derivatives of loop integrals. In these cases it is possible to proof gauge 
invariance analytically. 

5.2. Photon/ Gluon radiation 

The IR divergences can be removed using soft bremsstrahlung or by 
adding the corresponding 1 to 3 process with a massless particle (hard brems- 
strahlung). Soft radiation is proportional to the tree-level width but depen- 
dent on the energy cut /S.E of the radiated-off particle. It is automatically 
included in FC, see the formulae in 2^ . For an 1 to 3 process with a massless 



particle the three-body phase space can still be integrated out analytically. 
We have implemented this radiation by using self-derived generic formulae 
for all four graphs in Fig. [T] where every charged line can radiate off a photon 
(or a gluon for colored particles). The IR convergent total width is then 
given by 

r*°*^' = r(0o ^ P1P2) + r(0o ^ pm i/g) ■ 

For the simplest case, the decay into two scalars, r(0 — )• 0i02 l/g) is propor- 
tional to 



- 4|5(ti.eeP X [ g^mlloo + 9o9i{rnl + ml-ml)Iio + gfrnjlu 

I ^2 2|2\7-|22ri /2 2 2\r 

+ gog2{mQ -m^+ m^jho + 921^2^2 + 9i92[mQ - - m^Jhi 
+ 9oi9i + 92)Io + 9ii9o - 92)h + 92(90 - 9i)h] ■ 

The 'bremsstrahlung integrals' I are given in [26]. The integrals lij depend 
on log A, here A is the auxiliary mass for j/g. For the cases scalar — )■ fermion 
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+ fermion with one fermion mass zero (e.g. — )■ r+z/^), we have derived 
special formulae for the bremsstrahlung integrals. The other formulae can 
be found explicitly in the program code in the file bremsstrahlung. F. 

6. Program manual 

6. 1 . Requirements 

• Fortran 77 (g77, ifort77) 

• C compiler (e.g. gcc) 

• Loop Tools jiij 

6.2. About version 1.0 

• The CKM matrix is set diagonal 

• Real SUSY input parameters 

6. 3. Installation 

1. Download the file hfold.tar at 

http:/ /www.hephy.at/tools| 

2. expand the file, go to the folder hfold/SLHALib-2.2 and type 

. /configure 

make 

3. to create the Fortran code for hf old, go back to the folder hfold and 
type 

. /configure 
make 

4. To run HFOLD type 

hfold 
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6.4- The input file hf old. in 

1. name of the spectrum (SLHA format) 

2. Higgs boson = 1,2,3,4,5 

1 = /iO, 2 = //o, 3 = AO, A = H+, 5 = All 

3. contribution = 0,1,2 

= tree-level calculation 

1 = full one-loop calculation 

2 = SQCD (only diagrams with a gluon/gluino are taken into account) 

4. bremsstrahlung = 0,1,2 

= off, 1 = hard bremsstrahlung, 2 = soft bremsstrahlung 

5. resummation of bottom yukawa coupling = 0,1 

= off , 1 = on 

6. esoftmax 

cut on the soft photon (gluon) energy, if soft strahlung is used 

7. name of output-file 



7. Comparison HFOLD with HDECAY 3.53 and FEYNHIGGS 
2.7.4 

SPSla' point: 

In the following we show some results for the mSUGRA point proposed in the 



SPA project |25|, (Mi/2,Mo,v4o) = (250,70,-300) GeV, sign(/i) = +1, and 
tan P = 10. Our comparison with other programs is based on the same input 
file with the MSSM spectrum given in SUSY Les Houches accord from 24 



created by SPhenoS . Obeta. A list of available decay programs is given at 



http://home.fnal.gOv/~skands/siha/f In the following tables the Higgs 



bosons partial and total decay widths are compared to HDECAY 3.53 and 
FeynHiggs2.7.4. In FeynHiggs 2.7.4 the decays into fermions are at full one- 
loop level. HDECAY 3.53 [35|] has implemented higher order QCD and some 
EW corrections. Most of these corrections are mapped into running masses. 
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/i° HF-tree HF-SQCD HF-full FH 2.7.4 HD 3.53 

r*"*"' 1.9 3.0 2.8 3.2 3.7 



Table 1: Comparison of the total decay widths of the CP-even Higgs boson (in MeV) 

H° HF-tree HF-SQCD HF-full FH 2.7.4 HD 3.53 

Ytotai 0.8389 1.0171 1.0274 0.9890 1.0495 

Table 2: Comparison of the total decay widths of the CP-even Higgs boson 

A° HF-tree HF-SQCD HF-full FH 2.7.4 HD 3.53 

^totai ;l_2471 1.4405 1.5256 1.4183 1.4139 

Table 3: Comparison of the total decay widths of the CP-odd Higgs boson A° 

H+ HF-tree HF-SQCD HF-full FH 2.7.4 HDECAY 

Ytotai 0.7534 0.9057 0.8948 0.7875 0.9671 

Table 4: Comparison of the total decay widths of the charged Higgs boson i/+ 

BR-tree HF-tree HF-SQCD HF-full FH 2.7.4 HD 3.53 



bb 


0.8044 


0.0015 


0.0026 


0.0024 


0.0025 


0.0029 


TT 


0.1544 


0.0003 


0.0003 


0.0003 


0.0003 


0.0003 


CC 


0.0403 


0.0001 


0.0001 


0.0001 


0.0001 


0.0001 


Table 5: Comparison of the partial decay widths of the CP-even Higgs boson 




BR-tree 


HF-tree 


HF-SQCD HF-full 


FH 2.7.4 


HD 3.53 


66 


0.5546 


0.4652 


0.6262 


0.6216 


0.6283 


0.6466 


TT 


0.1058 


0.0887 


0.0887 


0.0914 


0.0983 


0.0909 


tt 


0.0549 


0.0460 


0.0631 


0.0564 


0.0607 


0.0937 


/VI A2 


0.0539 


0.0452 


0.0452 


0.0465 


0.0429 


0.0442 




0.0515 


0.0432 


0.0432 


0.0527 


0.0528 


0.0568 


Tlfi 


0.0212 


0.0177 


0.0177 


0.0184 


0.0183 


0.0095 


fif2 


0.0206 


0.0173 


0.0173 


0.0191 


0.0183 


0.0262 


-0-0 
A2A2 


0.0205 


0.0172 


0.0172 


0.0206 


0.0210 


0.0225 




0.0172 


0.0144 


0.0144 


0.0140 


0.0122 


0.0127 



Table 6: Comparison of the partial decay widths of the CP-even Higgs boson H' 
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BR-tree 


HF-tree 


HF-SQCD 


HF-full 


FH 2.7.4 


HD 3.53 


66 


0.3741 


0.4665 


0.6282 


0.6250 


0.6269 


0.6439 




0.1800 


0.2245 


0.2245 


0.2862 


0.2395 


0.2389 


tt 


0.0862 


0.1074 


0.1389 


0.1289 


0.1881 


0.1815 


-0-0 
A.xA.2 


0.0755 


0.0942 


0.0942 


0.0972 


0.0890 


0.0871 


A,2A,2 


0.0729 


0.0909 


0.0909 


0.1166 


0.0975 


0.0955 


TT 


0.0713 


0.0889 


0.0889 


0.0919 


0.0980 


0.0911 


T1T2 


0.0225 


0.0280 


0.0280 


0.0297 


0.0292 


0.0272 




0.0170 


0.0212 


0.0212 


0.0205 


0.0181 


0.0183 



Table 7: Partial decay widths of A° 



H+ 


BR-tree 


HF-tree 


HF-sqcd 


tb 


0.6171 


0.4649 


0.6170 




0.1712 


0.1290 


0.1290 




0.1203 


0.0906 


0.0906 




0.0809 


0.0610 


0.0610 



Table 8: Comparison of the decay widths 



HF-full FH 2.7.4 HDECAY 3.53 

0.5989 0.5060 0.6850 
0.1306 0.1194 0.1228 
0.0944 0.0922 0.0927 
0.0643 0.0630 0.0581 

the charged Higgs boson if+ 
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The screen output when running HFOLD is as following: 



J __ \ 



\ 



l__l 



l__ 



LI LLI 



\ /I. 



./ 1.0 



Higgs Full One Loop Decays by W. Frisch, 
H. Eberl, H. Hlucha 



error 

abort 173248520 
nslhadata 5504 



FF 2.0, a package to evaluate one-loop integrals 
written by G. J. van Oldenborgh, NIKHEF-H, Amsterdam 



for the algorithms used see preprint NIKHEF-H 89/17, 
'New Algorithms for One-loop Integrals', by G.J. van 
Oldenborgh and J. A.M. Vermaseren, published in 
Zeitschrift fuer Physik 046(1990)425. 



ffxdbO: IR divergent BO', using cutoff 1. 
ffxcOi: infra-red divergent threepoint function, 
working with a cutoff 1. 
Flags : 



Susyqcd calculation 



resuiranation of bottom yukawa coupling off 



Using onshell Higgs masses from : SPSlaprime . spc 



hard bremsstralilung on 
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the output SLHA will be written to : output. slha 



Decay- 


Table : 




Total 


width : 0.905677243 


H+ -> 


niu+ nu_mu 





H+ -> 


tau+ nu_tau 





H+ -> 


c sb 





H+ -> 


t bb 





H+ -> 


~chi_10 ~chi_l+ 





H+ -> 


~chi_30 ~chi_l+ 





H+ -> 


h W+ 





H+ -> 


~nu_eL ~el_2+ 





H+ -> 


~nu_muL ~mu_l+ 





H+ -> 


~nu_muL ~mu_2+ 





H+ -> 


~nu_tauL ~tau_l+ 





H+ -> 


~nu_tauL ~tau_2+ 






320441E-003 


/ 


BR 





00 


906196E-001 


/ 


BR 





10 


349004E-003 


/ 


BR 





00 


617035E+000 


/ 


BR 





68 


128997E+000 


/ 


BR 





14 


699102E-003 


/ 


BR 





00 


277014E-002 


/ 


BR 





00 


126903E-002 


/ 


BR 





00 


237453E-003 


/ 


BR 





00 


124702E-002 


/ 


BR 





00 


609862E-001 


/ 


BR 





07 


114467E-002 


/ 


BR 





00 
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